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Abstract 

In the framework of the chiral perturbation theory (xPT) we computed two- and three-loop 
corrections to the pion parton distributions which posses (5-function singularities at xsj — > 0. From 
this calculation one can conclude that in the region of small xsj ~ m% / {'^■k Ft^)'^ standard xPT needs 
QQ \ in the resummation of all orders in order to get correct description of the PDFs at small- a; region. 

■ We demonstrate an example of such resummation. 
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1 Introduction 



j dx q(x) = 1' y" ^ Qi^) = ^^2 ■ (2) 



Recently chiral perturbation theory (xPT) has been applied to parton distribution functions (PDFs) 
[U [21 m [3] in order to determined the dependence of these quantities on the mass of Goldstone bosons 
(rn^) and/or external small momenta and t in the case of generalized parton distributions (GPDs)). 
In contrast to standard xPT, in this case one deals with hadronic matrix elements of non-local light-cone 
operators (n^ — 0): 

0(A)=q-(iAn)7+{l, ir'^}q(-iAn). (1) 

1^ I Therefore apart from chiral expansion parameters (tott) small external momenta) one has the additional 

■ scale A - the distance characterizing the non-locality of the operator. Interesting theoretical question is 
to study the interplay of this additional scale with the chiral expansion parameters. In the present paper 

' we investigate this problem for the quark parton distributions in the pion. 

^ I For the case of the quark operators there exist two independent distributions corresponding to the 

QP . two possible values of the isospin I: I — defines the singlet distribution Q{x) and I — 1 non-singlet one 

I q{x), which are normalized as: 

(N 

l> 
O 

» ■ o ■ 

. Here is the momentum fraction carried by quarks. For convenience, the explicit definitions of the 

PDFs in terms of the matrix elements of the light-cone operators are presented in the Appendix. 
The standard chiral expansion of these PDFs implies: 

X 

Q{x) ^Q^°^{x) + a^Q(^\x) + alQ(^'>{x) + ..., (3) 
q{x) ^q<^"\x) + a^q<-^\x) + alq<-^Hx) + (4) 

where ~ (mTr/dTri^jr )^ {Ftt ~ 93 MeV is the pion decay constant) is the chiral expansion parameter. 
Important observation is that in the above formulae one assumes that the momentum fraction x has no 
chiral power counting, i.e. 

X ~ 0{al). (5) 
The one-loop non-analytic contributions have been computed in many papers [TJ [51 [31 11] . They read: 

Q'-^\x) = X ln[l/aj + 0(0^), (6) 
gW(x) = {g(")(a:)-J(a:)}ln[l/aJ + 0(a^). (7) 

We observe that isovector PDF q{x) has specific 5-singular behavior at small momentum fraction. This 
can be understood as following: for some small values of the momentum fractions x ~ the next-to- 
leading singular term formally is of the same order as the leading term a-^S(x) ^ O(a^). Moreover, higher 
orders of xPT may posses even more singular structures, like derivatives of the (5— function: S^"'^ (x). This 



1 



may indicate that in this situation the usual chiral expansion is not vahd anymore. The question is: 
How to write correctly the chiral expansion for x ^ { equivalently for the large light-cone distance 
A ~ l/ox) chiral expansion for this case controlled by the finite order of the standard xPT? 

It is clear that in order to construct the improved chiral expansion for the region x ^ one has to 
perform resummation of all (5-singular terms. Then such reordering of the xPT expansion may lead to 
generation of the some finite size function instead of singular 5{x) and therefore provides more correct 
chiral description of the small— x behavior. 

Let us formulate the above suggestion more accurately. We shall assume, that performing the chiral 
limit at fixed value of the momentum fraction x one obtains the standard chiral expansions ([3|) , ([4]) which 
can be rewritten as a sum of the regular /'"'^^(x) and (5— singular parts 

Q{x) = g-9(x) + S'{x)Di + S'"{x)D, + ... = Q-9(x) + E,>i, odd A S(^\x) , (8) 
qix) = r7-9(x) + S{x)Do + 5"{x)D^ + ... = q-^3{x) + E,>o,cvc„ A 5i^{x) , (9) 

where the regular contributions can be represented in the framework of xPT by series similar to p|l. 
(HI) and where each term in the expansions is represented by some smooth function of the momentum 
fraction x. The coefficients Di in front of derivatives of the (5— functions are some constants which also 
can be expanded with respect to small parameters a^^: 

D.^Df^+a,Df>+o^f^ + ...^ (10) 

We shall assume that PDFs are smooth functions in the chiral limit, i.e., for all values z = 0, 1, ... 

Df)=0. (11) 

Then the singular terms can occur only from the corrections induced by the loop diagrams of the chiral 
perturbation theory. The power of the derivatives is dictated by the symmetry properties of the given 
PDF. Recall that singlet (non-singlet) pion PDFs can be understood as antisymmetrical (symmetrical) 
function with respect to exchange x <-> —x being extended to the whole region — 1 < a; < 1. 

In order to check the assumption ^ or ^ it would be interesting to see that the singular contributions 
~ 5^^\x) can occur at least in the first few orders of xPT expansion. In next section we compute such 
singular contributions to the PDFs and find the explicit contributions of the type [a^ \a.{l / a^)]^'^^ 5^"^^{x) 
for the cases n = 1,2. We shall see that the mechanism of the generation of (5— functions is quite general 
and should work also at higher orders. In the section Discussion we briefiy discuss the main consequences 
which follow from representations ([8]) and ([9]) and possibilities to perform the resummation of the singular 
contributions. On the concrete example we demonstrate that resummed PDF is given by finite size 
function f{x/a^) which can be represented by the infinite sum of the (5— functions being expanded with 
respect to a^. 



2 The singular (^—function contributions in ChPT up to three 
loops 

Our goal is to investigate the higher orders of the chiral expansion in order to establish the existence 
(5— singular terms. We shall accept following strategy: we restrict our consideration only to leading 
non-analytic contributions. Such approach allows us to perform all calculations in the framework of the 
leading order chiral Lagrangian 

£(2) = i tr {d^Ud^U^ + x^U + xU^) , (12) 

for the pion dynamics and leading order chiral representation for the light-cone operators ([l]). For 
convenience, all necessary technical details and definitions are presented in the Appendix. The accuracy 
of our approach, despite we are working at higher orders, allows us to neglect the renormalization of the 
pion mass and chiral couplings, hence below we always assume their physical values. 
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Figure 1: One-loop diagram which generates 5— contribution to the pion PDF. The large blob denotes two- 
pion non-local operator vertex generated by corresponding chiral Lagrangian for the light-cone operator 




The calculation of the contributions of double chiral logarithms is well-known subject and have been 
already discussed in the literature, see for instance [5]. The main observation is that such terms occur 
due to the t/y— divergencies of the xPT diagrams and can be effectively computed if one uses the basic 
property of locality of the t/]/— counterterms. 

We shall use dimensional regularization with space-time dimension d = 4 — 2e in order to compute the 
bare matrix elements. For the one-loop diagram the coefficient in front of chiral logarithm is determined 
completely by the UV—pole, therefore it is enough to compute only the 1/e-pole part of the graph. This 
is clear from the general structure for the unrcnormalized matrix element: 

£e-^+-^ {p\OiX)\p) Inoop ^ {O^'^) = (^^) + ... = (13) 

(14) 

where the chiral scale fi^ ^ AttF.^. The interesting for us 5— term is generated from the diagram in FiglT] 
Let us describe shortly corresponding calculations. The expression for the diagram is given 

Gi = ^e[abc] a, ,(°)(/3) * / dk M(^P+ - [„4(fcp) + ...j ^ (15) 

where the dots in the brackets denote irrelevant (regular) contributions from the four-pion vertex and 
by the asterisk we denote the convolution integral with respect to (3. We also neglected various terms 
with e which can provide only the finite terms and denote dk = d'^k/n'^^^. Obviously, the integral in 
is quadratically divergent. Moreover, from the structure of the denominator and numerator one observe 
that fc+ in the argument of the 5— function can be ignored due to the rotation invariance and one obtains 

Gi = -ie[abc] S{x) / dk —— -3, (16) 

J ml [fc^ — ml\ 

where we used the normalization condition ([2]). The integral over the momentum k can be easily computed 
and the result defines the constant D^'' which reproduces the contribution with the 5— function from ([7]). 

The next important observation is very simple. One can expect that at the higher orders the non- 
trivial contributions to the coefficients Di can be generated in similar way if the four-pion subgraphs 
generates the scalar products (fc ■ pY with the highest possible power i in the numerator. Then the 
i -\- 1-loop graph for the forward matrix element includes the following master integral 

. u« . '--f!!^- -if-' '"■':!!:" + (17) 



^ J [k^-mV^ 



m-n 



This integral can be computed in the same way as before, but now one has to perform the expansion of the 
(5— function with respect to k+ in in order to obtain factor which is necessary in order to contract mo- 
mentum k in the scalar products. Hence the master integral can provide contribution with '-^ (5'-*'(a;)a^"*"^. 
At the same time, the integral is [/T^— divergent and therefore the result contains the contribution with 
the maximal power of the logarithm ^ a^-^^ ln*'''^[l/a^]. This is exactly the contribution which we want 
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Figure 2: Two- loop diagrams with the two-pion operator vertex contributing to the chiral expansion of 
the pion parton distributions. 



Figure 3: Two-loop diagram with the true four-pion operator vertex (without tadpoles) contributing 
to the chiral expansion of the pion parton distributions. The singular (5— functions are absent in the 
contribution from this graph, see explanation in the text. 



to compute. All this is an illustration of the main idea how the higher derivatives of the (5— function may 
appear at the higher orders. In order to perform the explicit calculation of the model independent loga- 
rithmic contributions one has to take into account the presence of the subdivergencies in the multiloop 
graphs. The idea of the method is quite general and based on the locality of the f/F— counterterms. For 
the two-loop case it was discussed in Ref. [S| but can be easily extended to arbitrary order of xPT. For 
convenience, we provide below the technical analysis for the two-loop calculation of chiral logarithms in 
our case. 

The renormalized matrix element of an operator to two-loop accuracy can be written as: 



e 



(18) 



Let us clarify the rhs of eq. (HH]). The first term (O'^^^) is the bare matrix element (two-loop graph). The 
second term i (O^^^) denotes subtractions of the one-loop subdivergencies associated with the operator 
vertex or pion Green functions. The one- loop reduced matrix element (O*-^^) includes effective vertex 
from the corresponding one-loop counterterm. And the third term is proportional to the tree level matrix 
element (0(°)) of some effective operator. This term corresponds to the local counterterm, which removes 

the total divergency. The dimensionless combination t = is introduced for convenience. To two-loop 
accuracy the structure of divergencies of each contribution in eq.((T8]) has the following form: 



(19) 



= ^.92 + -gi + .90, {0^'^) = -hi + 



where functions gt and hi are polynomials of external momenta and the pion mass. Inserting these 
expansions into (fT8|) one obtains 



1 _ , . ._2, 2 . 1 



{0)j^ = ^g2[l + 2e \nt + 2e^lnt] + — hi 



1 + e lnt+ -^e'^^^^t 



\[g2+ hi] + -[2g2 + hi]\nt+ -[...] +\n^t 



2,92 + 2 ^1 



(20) 
(21) 



local terms 



Now, the non-local 1 /e divergency proportional to In t must cancel (locality of counterterms) , that pro- 
vides the relation: 



252 + ^1 = 0, ^ 52 = hi. 



(22) 
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Figure 4: One-loop sub-divergencies of the two-loop diagram Gi and corresponding counterterm. The 
filled circle is an operator counter, filled square is the 4-pion Green function counterterm. 



i.e., the residue of the is defined in terms of the one-loop diagrams with insertions of one-loop 

counterterms. This defines completely the coefficient in front of double logarithmic contributions: 

2g2 + hi = ~hi (23) 

Therefore, in order to define the contribution of the type ^ In^(mJ) we can compute only pole in 
the diagrams with insertion of the one-loop counterterms. In addition, we have to select such structures, 
which can produce (5— functions. From the observation considered above, one can expect that the two-loop 
diagrams can contain master integrals eq. (|17p with the power i = I, hence such diagrams are relevant for 
the computation of the coefficient Di in front of S'{x) from ([8]) (7 = 0). 

At the NNLO of xPT we have set of the two- loop diagrams presented in Fig. [2] and Fig. \3\ Diagrams 
G3 — Gg are insertions of the tadpole loops to the one-loop diagrams. From simple analysis one can 
conclude that these diagrams do not contain the master integral like (fT7|) with the maximal power i = 2 
for scalar product (fc • p) because of dimensional reasons: the tadpole integral is proportional to the pion 
mass ~ and therefore reduces the possible power i by one unit. It means that these diagrams have the 
same singular structures as the one-loop graphs, i.e., they can have only S{x) terms without derivatives. 
The more singular terms may appear only from the first two diagrams Gi 2- 

Consider these graphs. There are two one-loop divergent subgraphs in these diagrams: operator 
vertex subgraph and pure pion subgraph. They are depicted in Fig. [H There are two corresponding 
counterterms: 

(OW) = it- + (0«), (24) 

Computation of these contributions provides following result 



,2 



(25) 



No new singular structures of the isovector PDF q{x) appear, as the isovector PDF is symmetric. 

Let us briefly discuss the diagram on Fig[31 This diagram includes the four-pion operator vertesfl- 
Such vertex can be generated either from the operators (|17|) and (|48p or from the multiparticle terms 
which we do not write explicitly in Eqs. (|47l48p . These additional multiparticle terms of the chiral 
operator include new chiral constants and do not generate two pion vertices. Their expansion begins 
from the four-pion contributions as a minimum. We checked that the diagram Fig. [3] does not produce 
the contributions with the (5— functions in the case when the vertex is generated from the operators (|47|) 
and This can be seen from the fact that one pion from the operator vertex describes the external 

state and therefore provides the non-trivial argument for the distribution of the fraction x that excludes 
occurrence S{x). The contributions of the multiparticle terms with the new chiral constants have been 
ignored for the same reason. 

Now let us briefly discuss the general mechanism of generation (5— singular terms. It turned out, that 
such terms appear only from the diagrams with four-pion counterterms and only in the diagram Gi. 
We checked that the operator vertex counterterm, provides only simple S{x) without derivatives. The 



^We do not consider here the operator vertices with the tadpole loops as the true multiparticle vertices 
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Figure 5: Three- loop diagrams with the four-pion operator vertices 




Figure 6: Three-loop diagrams which may produce (5— function contributions. 



term with the pion counterterm exactly reproduces the master integral ([T7|) with i = 2. It also naturally 
explains the absence of S'{x) at one loop: the four-pion vertex is quadratically divergent and in the 
corresponding master integral can occur only simple (5— function structure considered in eq. (|15[) . Note 
also that the 5— producing master integrals have always maximal value of the f/l^— divergency index, i.e., 
power divergent terms play the key role in the generation of the (5— singular contributions. 

Previous analysis can be easily extended to the three-loop diagrams. Following the method discussed 
above one can again reduce the calculation of the a|ln^[l/a;^] contribution to the reduced one loop 
diagram with effective vertex corresponding to the insertion of the two-loop counterterms. 

We again have the subset of the graphs with the four-pion operator vertices shown in FigO In the 
case of isospin 1=1 such vertices appear only from the multiparticle part of the effective twist-2 operator 
with the new chiral constants. Therefore such contribution without loss of generality can be considered 
separately. In this paper we shall ignore them what, however, cannot change our main conclusions. 

The set of the "two-pion" relevant graphs is given in Fig. [HI In Fig. [7] we show the two-loop subgraphs 
for the case of diagram Gi and resulting one-loop graphs with effective vertices. We do not show all 
diagrams with tadpole loops and wave function rcnormalization. It is clear that such contributions can 
produce only 5— functions which occur in the one- and two-loop graphs. But now we are interested in 
the new, more singular terms with the second derivative 026^^'' (x) in Eq. The main steps of the 
calculation are the same as for the two-loop case. Again the singular contributions occur only from the 
reduced diagrams with effective vertex generated by the four-pion subgraphs ( the upper line in Fig. [7]). 
The substitution of the effective counterterm vertex reduce the calculation to the master integral (|17p . 
Only the diagrams Gi, G4 and G5 provides non-trivial contribution to the coefBcient D2- The three-loops 
result is 

D2=alln'[l/a,] (x^) (- ^) . (26) 
where (x^) denotes the moments of the PDFs in the chiral limit 

{x"^^) = [ dxx"''\^°'>{x) (27) 
Jo 

It is obvious, that given analysis can be continued further order by order in xPT. We do not find 
any special reasons that could suppress contributions with the higher derivatives S'-'^\x) in ([8]) and ([9]). 
Hence, we may conclude that structure of the chiral expansions for pion PDFs proposed in formulas ([5]) 
and ^ is established. 
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Figure 7: Diagram Gi (left) and two reduced graphs with effective vertices (right) which have to be 
computed in order to define (5— function logarithmic contribution. By the dashed hne we show the two- 
loop subgraphs which has to be substituted by its J/V— counterterm. 



3 Calculation of the singular terms in the pion GPDs 

In this section we discuss shortly the results for the singular contributions to the pion GPDs, counterparts 
of the (5— contributions to the pion PDFs, see eq. ([37|) . We shall consider, for simplicity, the limit 



i 0, C 7^ 0. (28) 
Logarithmic one- loop calculation gives [21 |4j : 



H'=\x, Ok^o, t^o = H'-'{x, ( 1 + «x In [l/«x]) - «x In [l/«x] 'f [^j ■ (29) 

Here the second term in the forward limit 0) tends to the S{x) contribution for the pion PDF 

«(H<0. /xv^„^j^j (30) 



In order to obtain above result we used that pion distribution amplitude ipTriz) in the chiral limit, due 
to the soft-pion theorem [5], can be represented in terms of double distribution as 

ipAz) = F'=\p,a) *Sia + f3-z), (31) 

The isoscalar pion GPD does not have logarithmic contribution at the one-loop order. Calculations 
of the singular higher loop contribution goes along the line described in previous section, the final result 
is as follows: 



H^-'ix,0 = <t(x,0 + I 4ln^ [1/aJ d(^^], (32) 



Here D{z) is the so-called D-term [7 , which in the chiral limit, due to the soft-pion theorem [6 , can be 
expressed as: 

D{z) = 2F^=°{P,a) 'f'6{a + P-z), (34) 
and function i^{z) is expressed in terms of double distribution as follows: 

(z) = lim —r,F^=\/3, a) * p 5 (a + p7] - z) . (35) 

One can easily check that in the forward limit (^ 0) the equation ([5^ is reduced to the result for the 
pion PDFs (i?7| . 

The singular terms of GPDs ([5^ result in the contributions to the amplitudes for hard exclusive 
processes of the following type: 
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71=1,2,3... (36) 



This clearly demonstrates that any finite order of chiral corrections to physical observables at ^ ~ 
a;^ln(l/a^) gives contribution which is not suppressed by small chiral parameter. Therefore in this re- 
gion of kinematical variable ^ one needs to perform infinite order resummation of the chiral perturbation 
theory. 



4 Discussion 

In present paper we have demonstrated that chiral expansion of the pion PDFs includes the contributions 
of the 5— function and its derivatives. We have also computed the leading non-analytic contributions to 
the coefficients Di 2,3 in ^ and (O. The results read: 



Q{x) = Q-'^a^x) + ( -- ) al In^ [1/aJ 5'{x) -f ©(a^ ), (37) 



q{x) = g-^x) - In [1/aJ 5{x) + {x') W [1/aJ 5" {x) + 0{a\), 



(38) 



It is also clear that higher order terms of the chiral Lagrangian for pion dynamics /^6,8... a-nd for the 
light-cone operators can also produce (5— singular terms but they arc suppressed by powers of the small 
chiral parameter with respect to leading chiral logarithms (~ [a^ ln(l/a;^)]) by additional powers of a^. 

From the obtained results we can conclude that in the region of small momentum fraction x application 
of the standard chiral expansion faces with the difficulties. In this region one cannot truncate the 
perturbation expansion and hence to perform some sort of resummation. This reorganization of the 
chiral expansion may combine the singular contributions to some generating function f{z): 

Y,Dne^5^''\x) = fix /e) (39) 

where small parameter e — a^ln[l/a^]. The function f{z) is some stable in the chiral limit function, 
which, probably, can be defined in the whole real axes —00 < z < 00 and generates contributions to the 
moments in the form 

1 POO 

dx f{x /e) « £" / dz z"-i f{z) ^ £" /„ = £" (-l)"(n - 1)! A^-i- (40) 

1 J -00 

The obvious conclusion which follows from this discussion is that any finite order calculation of the PDF 
moments can not provide enough information in order to reconstruct PDF in the whole region of variable 
X. One has to perform resummation of singular (5— function terms but this is equivalent to calculation of 
all orders xPT. 

Of course, such resummation is very difficult task. Even if we restrict our consideration only by the 
leading non-analytic terms we must compute straightforwardly order by order all terms which is not 
possible. Hence one has to introduce some simplifications or models in order to achieve the goal. As an 
example, let us provide here the result of such simplified resummation in the leading order of the 1/A^ 
expansion, where N is number of pions {N = 3 for the real two flavor QCD). This expansion is obtained 
by the generafization of the SU{2) x SU{2) ^ 0(4) chiral model to the 0{N + 1) model. The latter 
model can be solved in the large N limit. The result of the calculations for the isovector PDF isH: 

9 dz 
q{x)^qr,g{x)~4- 0{Ne-\x\) - q^"Hz). (41) 



N 



x\/{Ne) ^ 



[Note, that parametrically e = l/2a^ In [1/a^] ~ ^/N]. This example of chiral resummation shows clearly 
that calculation of any finite order of xPT contributions to the Mellin moments of PDF docs not allow to 



^Details of the calculations and discussion will be given elsewhere [8] 
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restore correctly the distribution function for x < 3/2a^ In [l/a,^] (numerically for the physical pion mass 
X < 0.09 that is not very small). For the chiral extrapolation of the lattice data for PDFs this implies 
that usage of the moments without resummation (see e.g. [9]) can provide us the incomplete information 
about parton distributions in the region of a; ln[l/ap^]. In other words, one can use xPT in order to 
compute systematically the local moments of the PDFs but one cannot use these moments in order to 
describe the PDFs in the region of the small momentum fraction where all orders of the xPT contribute 
to the same accuracy. 
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6 Appendix. Twist-2 operators and the matrix elements 

In this section we briefly describe the definitions and some technical details used in the paper. We 
introduce two light-like vectors n, n: 



n =n — 0,n ■ n ~ 1, a-f- ~ a ■ n. 
There exist two QCD quark light-cone operators of twist-2: 

PB. = lil-To), Pl = ^(1 + 75), 



(42) 



(43) 



[OR]f, = qg (^A n) ^+Pn qf (-lA n) , (44) 
[OL]fg = % (iA n) qf (-lA n) , (45) 

where indexes /, g stand for flavor. These operators transform under global chiral rotations as 

Ol ^ VlOlVI , Or ^ VrOrVI . (46) 

In xPT these QCD operators are described by effective chiral operator with unknown chiral constants. 
In the pure pion sector {U as usually denotes pion field) one finds [21 SI 13] 



U ( ^An") n. a C/t f ^An 



[/t ( ^Xn\ n-du( ^Xn 



where by asterisk we denote the integral convolution with respect to /3 and a: 
T{P,a)*0{(3,a) = d(3 da T{f3, a) 0{f3, a) . 



fg 



fg 



(47) 
(48) 

(49) 



Here J- {(3, a) represents the real generating function of the tower of low-energy constants and df^ denotes 

combination of derivatives — d^. The non-local structure of the operators (|T7)) and (HS)) implies that 
the light-cone distance A has chiral counting A ~ 0{p~^) as suggested in ^2). Then the combination A9+ 
is dimensionless that provides the non-locality with respect to light-cone direction. Important to note 
that given Lagrangian is not complete: the expressions (|47|) and (j48|) describe correctly only operators 
vertices with two attached pions (including the tadpoles loops) . The low-energy constants .F(/3, a) are 
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characteristics of the structure of the pion, they are not determined in the effective field theory. According 
to the isospin 7 = 0,1 one can construct two independent functions: 

F^=°[/3,a] = [-A a] [/?,«]), (50) 

F^=^[P,a] ^ ^{T[-f3,a]+T[(3,a]), (51) 

which are convenient to describe the pion matrix elements. Pion PDFs are defined as 

^e-^P^-^{n\p)\tr[T^OL+Bm\^''{p)) ^ ^^e[abc]q{x), (52) 

^e-^P+^^{n\p)\tr[OL+RiX)]\n^ip)) ^25-'Q{x) (53) 
that in the chiral limit can be written as 

r-l-|/3| 

da F'=°{(3, a) = [(?(/?) g(/3) - 0(-/3) g-(-/3)] = Q{(3), (54) 

-1 + 1/31 
1-1/31 

daF'='{(3,a) = 0(/3) q(/3) + 0(-/3) g(-/3) = <z(/3). (55) 

-1 + 1/31 



SumR The more general functions GPDs are defined as 



dX 
2^' 



(7r^(p')|tr[r=Oi+fl(A)]K'^(p)) ^ Aie[abc]Hi=\x,S„t), (56) 



^e-^^+-^(^^(p')|tr[OL+fl(A)]K'^(p)) ^25^^H'=^{x,i,t) (57) 

with P =\{p + p'), ^ = ~WW)T' ^ ^ (P' ~ Py^- I" forward limit ^ ^ 0, t ^ 0: 

H^=\x,0,0) = q{x), H^=°{x,0,0) ^ Q{x). (58) 
And the pion distribution amplitude is given 

6-'F^ ,p^{u) = I / |Ae-»(p-«)V2 (^a(p) ^^bQ^^^^X)] I 0) , (59) 
Defined distribution amplitude satisfies normalization 

j du 4>^{u) = 1. (60) 
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